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Abstract 



We compute the low energy threshold corrections to neutrino masses and mixing 
in the Standard Model (SM) and its minimal supersymmetric version, using the 
effective theory technique. We demonstrate that they stabilize the renormalization 
group (RG) running with respect to the choice of the scale to which the RG equation 
is integrated. (This confirms the correctness of the recent re-derivation of the RGE 
for the SM in |hep-ph /0108005 . ) The explicit formulae for the low energy threshold 
corrections corrections can be applied to specific models of neutrino masses and 
mixing. 



1 Introduction 



There is at present a strong experimental evidence for neutrino oscillations. Their most 
natural explanation is the existence of neutrino masses. Neutrino masses can be incorpo- 
rated in the Standard Model (SM) or its supersymmetric extension (MSSM) by adding 
to the Lagrangian the non-renormalizable dimension-5 operator [|I| 

A/:sM = -^C""^ {tMt) {ei.mf) + H.c. (1.1) 

in which A, B = 1,2,3 label generations, if = {u^, e^) are the Weyl spinors transforming 
as doublets of SUl{2), is the Higgs doublet with hypercharge +1/2 and €21 = —612 = 1, 
en = 0. After the electroweak symmetry breaking (|1.1|) gives the neutrino mass matrix in 
the form 



(where v is the vacuum expectation value of the neutral component of the Higgs field 
which is diagonalized by the unitary rotation z/^ U^"'Ua. The elements of the matrix U 
determine the neutrino oscillation probabilities and are, therefore, probed in the neutrino 
experiments. 

The operator (|1.1| ) appears in the low energy effective theory as a result of integrating 
out fields of an underlying theory describing physics at some high energy scale. Thus, it 
is supposed to be generated at some high scale Mp (much higher than the electroweak 
scale Mz). Therefore obtaining reliable predictions for neutrino masses and mixing angles 
requires solving the renormalization group equation (RGE) for the Wilson coefficient C^^ 



d r^AB T^j^AB I ,^ r ,2 y^AB I r^AB„,2 1 /i q\ 

— C = KG +K\y^^C +C ?/g^J (1.3) 

where t = (l/lGvr^) \ia{Q/Mz) and y"^^ are the Yukawa couplings of the charged leptons.Q 
In the SM K = -3/2 and K = -3g^ + 2 Efcrmions N^^^yj where N^^^ = 3 for quarks and 1 

2 0^2 I „.2 



for leptons; in the MSSM k = +1 and K = -6gi - 2^^ + GEaV^ 



The solution to eq. (|1.3|) 



C(Q) = IkJC{Mf)J (1.4) 



^Without loss of generality throughout the paper we work in the basis in which the matrix of the 
charged lepton Yukawa couplings is diagonal; the Wilson coefficient C^^ is therefore assumed to be 
given in that basis too. The normalization of A is fixed by the Higgs self interaction: £scif — —^{H^H)^. 
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where J = diag(/e, I^, It) and 

rtQ 



IK 



exp (- ^ K{t')dt' 







/e, = exp yiit')dt']^l-ri (1.5) 



with tq = (l/lGvr^) In^Mp / Q) , gives C((5) at the electroweak scale Q ~ Mz in terms of 
C{Mp). The RGE ( |1.3| ) was analyzed in many papers to see how much the initial 
pattern of neutrino masses and mixing angles generated at the scale Mp is modified by 
quantum corrections involving large logarithms hi{Mp/Mz) ^1. In particular, it has 
been found ||^ that eq. ( |1.3| ) exhibits a nontrivial fixed point structure. In many interesting 
cases (e.g. for degenerate or partially degenerate neutrino mass spectrum) that structure 
leads to the pattern of mixing angles that is not compatible with the present experimental 
indications (i.e. bimaximal mixing and small f/13 matrix element). It has been however 
pointed out p, ^ |10] that in the MSSM the so-called low energy threshold corrections 
which were neglected in previous analyses P] can in some cases be more important than 
the RG evolution and can change qualitatively the pattern obtained by solving eq. ( |1.3| ). 

In this paper we compute these low energy threshold corrections both in the SM and in 
the MSSM. We first show that they stabilize the results of the RG running with respect to 
the choice of the low energy scale Q (to clarify the points raised in the recently published 
paper [11] we demonstrate this explicitly in a pedagogical way) and asses their magnitude 
and dependence on the parameters of the MSSM. 



2 Standard Model 

In this Section we calculate one-loop corrections to the neutrino mass matrix in the SM. 
Our starting point is the SM Lagrangian (see e.g. ref. |T^) supplemented with the 
non-renormalizable term ( |1.1|) . All parameters of this Lagrangian are understood to be 
running parameters renormalized at the scale Q ~ Mz- Also the Wilson coefficient C^^ 
of the dimension 5 AL = 2 operator ( |1 . 1| ) is a renormalized parameter of the effective 
theory Lagrangian. Integrating its RGE ( |1.3| ) from the high scale Mp down to some scale 
Q ^ Mz resums potentially large corrections involving ln{Mp/Q) to all orders of the 
perturbation expansion. However, since the low energy scale Q is not a priori determined 
by any physical requirement (apart from the condition Q ~ Mz), the neutrino masses 
and mixing angles computed in the tree-level approximation from the Wilson coefficient 
C^^{Q) do depend (albeit weakly) on the actual choice of Q. This dependence can be 
removed by computing masses and mixing angles in the one-loop approximation in the 
MS scheme with the same renormalization scale Q. 
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Since the neutrino masses are orders of magnitude smaller than the electroweak scale, 
the calculation of the low energy threshold corrections is technically most easily achieved 
in the effective theory approach. At the scale Q ~ Mz all gauge and Higgs bosons are 
integrated out and the effective theory valid below the electroweak scale is constructed. 
In this low energy theory the one-loop neutrino mass matrix (^mJ,^^°°P^ is given by the 



tree- level term ( |1.2|) of the SM plus the one- loop (threshold) correction Am,y. The latter, 
apart from having leading In Q dependence that exactly matches the In Q dependence of 
the tree level mass {m^J'^'^)^^ (Q) ( |1.2| ), can also contain nontrivial Q-independent pieces. 
Writing the SM Higgs doublet as 

we get from (|1 . 1| ) the neutrino mass term and various interactions (we write down only 
those which will be relevant for us): 

A£sM = -i(m~)'"'i'^i'B (2.2) 

In principle to get the Feynman rules for neutrino mass eigenstates Ua one has also to 
rotate the neutrino fields ua — ^ U^^'Ua (where U^"' is the matrix diagonalizing {m}J'^^)^^) . 
This step is however unnecessary for our purpose since we will use everywhere massless 
neutrino propagators on internal lines.0 We can therefore compute directly the corrections 
to the tree level mass matrix C^^ in the basis in which it is not necessarily diagonal. 




Figure 1: One-particle irreducible threshold corrections. 

The strategy is now to integrate out heavy fields: W"^, 0° (as well as Goldstone 

bosons and G^) and to construct the effective Lagrangian valid below the scale Q ~ 

^Taking non-zero neutrino masses into account in propagators would amount to including effects 
(where M is the mass scale of the heavy neutrino states). Since we do not consider operators of dimension 
higher than five resulting from the seesaw mechanism we cannot compute effects consistently. 
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Mz- Up to terms of order OijP' /M^) where p ~ nii, <^ Mz is the external four momentum, 
one-loop effects (shown schematically in fig. |I]) of the heavy fields present in the full SM 
have to be simulated by the corrections 5z^^ and 5m^^ in the effective theory Lagrangian 



+ ... (2.3) 



2 

Redefining the neutrino fields to get their kinetic term canonical and using bz^^ 
-S^^(O), bm^^ = ^i^{0) (where S^^(p2) and S^^(p2) are defined in fig. I) one gets 

A'B / , \AB' 



(Am,)^^ = If,^ (m'r) + (m'r) if'B (2-4) 



where we have split 

E;:,^(0) = 1a^'^ {mTf + 1 {mTf A^'^. (2.5) 
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Figure 2: Contributions of the charged Goldstone boson. Heavy dots indicate vertices 
arising from the operator ( |1.1D . Crosses indicate fermion propagators with a helicity flip 
(i.e. with the fermion mass in the numerator). 

Goldstone bosons G^ contribution to Sy^ and are shown in figs. |^ and b, c, 
respectively, where the dots denote interaction vertex originating from A£sm given in 
eq. (|2.2|) . In a general i?^-gauge and working in the MS scheme one finds 

1 

2' 



s;^^(o) 



4M 



C^^So(eA, G±) + C^V So(eB, G 



(2.6) 



We have used the abbreviated notation for the standard two-point function 



5o(l,2) = 5o(0,mi,m2) 
5^(1,2)= ^Bo{p\rm,m2] 



1 



p2=0 



-1 + 
1 
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1 1 
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1 1 
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= In — - 
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where mi and m2 are the masses of particles 1 and 2. 
The boson exchange contributes only to Sy: 
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- M^)B',ieA, W^) + 5o(eA, W^) + 1 



+ 



di.AB ^ 

4 Ww 



ml-Ml,)B'^{eA,W^)-B^{eA.W^) -2i?o(e^,lV^) (2.7) 
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Figure 3: Contributions of Z° and neutral scalars. 
Together, and contributions give 
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AB 92 



4 



5o(eB,W^^)- U 



m 



w 



Bo{eB, 



(2.8) 



In the limit of massless neutrinos on internal lines, the contribution of exchange 
to arising from the diagram shown in Fig. |^a can by obtained from (|2.8| ) by setting 
there m^^ = and replacing W^{G^) Z°(G°), {gl + gl)/2. This gives 



jZ°{i) _ xAb 91 + 9y 

^AB — ^ o 



l-M|i?;(z/,ZO)+ezi?o(^,G°) 



(2.9) 



(we do not write any index on v to stress that neutrino masses are set to zero in the Bq 
functions). The diagram shown in Fig. ^d, which arises due to the non-zero Majorana 
mass insertion, contributes to S^^(O). It gives 



.z0(2) _ .ab92+9y^ 

^AB 



i?o(z^,Z°) 



4M| 



a(G°) 



4M| 



a(Z° 



:2.10) 
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Finaly the exchange of and (p^ in the diagram shown in Fig. ^ gives 

1 



tAB 



6 



AB 



2^2 



(2.11) 



where IGvr^a 



■1 + ln(m^/(5^)] is another standard loop function. The dependent 



part of this contribution cancels the dependence of ( p.lOp . 




+ 




W^,Z^, ghosts 



+ 




fermions 



Figure 4: Tadpole diagrams. 



The combined contribution of ( |2.8[ - p.llD is still gauge dependent because one has to 
include the contribution of tadpole diagrams shown in fig. ^. They give 



'■AB 



-6 



AB 



Mlv 



^AB 



■^AB 



IGtt^ 



1 - In 
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1 - In 

2 
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Ml' 



It is easy to see that the and ^z dependence of ( p.8| ) and ( p.9| ) is canceled out by 
eq. (111). 



To check that [rri],' 



loop 



AB 



is independent of the renormalization scale Q we must 
recall the RGE for the vacuum expectation value f^. Since in this approach f ^ is merely 
an abbreviation for —2m?/\ where m? and A are the (negative) mass squared parameter 
and the self coupling of the Higgs doublet, respectively, we have 



d n 

dt 



1^2 Id/ 
— m — ——A 

, m'^ dt Xdt I 



(2.13) 



where IQn'^t = In Q. The RG equations for and A read: 

d o 



—m 
dt 



m 



9 2 2 

-^92 - ^9y 



6A 
9 



2T 
3 



|a = 12A^ - {9gl + gl)\ + \gt + \glgl + \g'y + 4AT 



4T2 
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where T = E/ N^y} and T2 = E/ N^yj. Combining eq. (^) with the RGE (^) for 
C^^ we have therefore: 



C^^'iQyiQ) 



AM 



4M 



V 



Q 



1 /9 4 3 2 2 



A V4 



34 .rr.\..^Q 



9l + -^gWy + 4^?y - j In ^ 



(2.14) 



It is now easy to see that the leading In Q dependence of the tree level neutrino mass 
matrix {vn}'^^^^^^ (Q) cancels out with the explicit InQ dependence of the 1-loop correction 
( [^.4| ) to it. In particular, the 1/A terms in eq. ( [2.14| ) cancel the InQ dependence of the 
W^, (and their ghosts) and fermionic tadpoles in ( p.l2|) . 

It is also possible to consider v not as the tree-level VEV of the Higgs field but instead 
as the minimum of the full 1-loop effective potential. In such an approach there are no 
tadpoles ^ but since the effective potential gives finite v only in the Landau gauge, = 

= 0, the contributions (|2.8| - |2.11| ) must be taken in this gauge too. In this approach 
the RGE for f i_ioop is no longer given by eq. ( |2.13| ) but instead is determined from the 
anomalous dimension (also taken in the Landau gauge) of the Higgs field operator: 



d 
dt 



_ 2 2,3 2 

1-loop — "^i-ioop {'^32 ' 



2T 



(2.15) 



AB 



It is then easy to check that again the explicit InQ dependence of (jnl obtained 
from eqs. ( |2.8| - |2.11| ) in the Landau gauge cancels against the In Q dependence of the tree 



level mass matrix {1/AM)C^^ {Q)vl_^^^p{Q). 

In practice, the difference between the two approaches (which formally is a higher 
order effect) is not seen when v"^ (or fi_ioop) is expressed in terms of the physical Z° 
boson mass. For example, in the first approach one has 



4(M2 



ZJph 



92+9y 



ZJph 



(2.16) 



where Ilzz{M'^,Q) is the 1-PI self energy of the Z° boson computed for = M| and 
renormalized in the MS scheme with the renormalization scale Q. 



■^Using V determined from the full 1-loop potential is equivalent to saying that one expands the sym- 
metric Lagrangian around some initially unspecified v and determines the value of v from the requirement 
that the tree level tadpole (arising from a term in the Lagrangian that is linear in the Higgs field) cancels 
the 1-loop one. 
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Neglecting terms of order 0{m^^/M^^) and higher the final formula reads (tadpoles 
have canceled out) 

rAB _ " J „,2 I 3 Myy \ 

4 ) 

2fizz(M|,g)l 





- »7r' 



{Ml 



^Z)ph ) 

where we have adopted C,w,z = 1 and v"^ in the tree-level neutrino mass matrix is now 
given by 

TTttEM 

where and uem are the sinus of the Weinberg angle and fine structure constant, respec- 
tively, in the MS scheme and at the renormalization scale Q (for which one can take M^).0 
Note also that the factor — |ln^^^ in the first line of ( |2.17| ) confirms the correctness of 
the recent re-derivation of the SM RGE. 

From our discussion it should be clear that putting a particular emphasis on better 
stability with the renormalization scale Q of the product v'^C^^ (or of its eigenvalues) as in 
ref. [0] makes no sense in the quantum field theory. The physical neutrino masses defined 



as the poles of the propagators (or, in the one-loop approximation, as the appropriate 
coefficients of the effective Lagrangian ( ^.31 ) ) do not depend on the renormalization scale 
Q. On the other hand, by themselves big changes of C^^ during the RG evolution 
between Mp and Mz do not signal any instability. They reflect only the importance of 
the resumation of large logarithmic contributions (y^ In^Mp/Mz))"^ where n = 1, 2, ... in 
order to get reliable results for neutrino masses in terms of the Lagrangian parameters 
deflned at the scale Mp. Finally, let us notice that v'^{Q) disappears from the flnal 
formula ( |2.17D for neutrino masses. Therefore, the question whether v"^ is considered as 



an abbreviation for — m^/A (in which case its variation with Q is very rapid) or as the 
minimum of the full one-loop effective potential is inessential for the stabilization of the 
results for physical neutrino masses. 



"'In the standard way (see eg. Q) s'^{Mz) and cxEMiMz) can be expressed in terms of measurable 
quantities: a^M measured in the Thomson scattering and by the Fermi constant Gp- 
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3 The MSSM 



In this section we calculate one-loop corrections to neutrino mass matrix in the MSSM. We 
will see that, apart from stabilizing the results obtained from the RG analysis with respect 
to small changes of the final scale Q, they contain also InQ independent terms which can 
be more important than the RG evolution. In some situations |TU|, Q they can change 
the pattern of mixing and lead to relations between the mixing angles different from the 
one obtained at the infrared fixed point of the RGE (|1.3| ) p. We will use the notation 



and conventions of ref. in which the Feynman rules resulting from the renormalizable 
part of the MSSM (i.e. without the higher dimension operators) are collected. 

The neutrino masses arise in supersymmetric models from one of the dimension 5 
operators obtained by adding to the Lagrangian supersymmetric non-renormalizable terms 

AC oc C^^ / d'9 {e.,Ht^Lf) (e^.i^/^^Lf ) + H.c. 

= -^C^' [{^^^Hthf) )/f ) + {e^.hT'Lf) {e.M^^ Lf) (3.1) 

+ terms with auxiliary fields (3-2) 

where capital letters with a hat denote superfields and capital letters and lower case letters 
denote their scalar and fermionic components, respectively. In the second line we have 
fixed the normalization so that the first term coincides with the operator (|1 . 1| ) in the SM. 

Expressing the initial fields in terms of the physical ones as in ref. []T^ one gets the 
following terms (we write down only those which are relevant for our calculation): 

f-iAB f^AB 

- jj^Z'^Z'^H^H^^Bi^A + -^Z]^Z'jHl,Hl,UBUA (3.3) 

V2M ^ ' 2V2M ^ + kAj n 

As previously we use Weyl spinors here. From eq. ( p.3|) the necessary additional Feynman 
rules can be easily obtained. 

The contribution to the quantity I^^ defined in eq. ( p.4|) of the and bosons 
is the same as in the SM and is given by eq. ( [^.8| ).p| The contribution of (arising from 

'"^Strictly speaking, in supersymmetry one has to use the DRED scheme instead of DIMREG used 
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diagrams similar to the ones shown in fig. ^ is 
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2 M, 



(3.4) 



w 



where tan/3 = Vu/vd is the usual ratio of the VEVs of the two Higgs doublet if*^"-* and 

The contribution of Z° is the same as in the SM and is given by eqs. ( |2.9| , p.l0| ) while 
the contribution of neutral scalars (p. 11 ) is in the MSSM replaced by 



rscalars 



sin^ a a{H^) + cos^ a a{h^) - cos^ /? a(A°) - sin^ p a{G°)] (3.5) 



Since v"^/ sin^ P = '^u'^'^d ^^^^ term of eq. ( p.5|) cancels the dependence in eq. (|2.10|) 
as in the SM. The dependence of ( |2.8|) and (|2.9| ) on and ^z, respectively is again 
canceled out by the tadpole diagrams with G"^ and loops.Q 
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Figure 5: Contributions of charginos/charged sleptons (a-c) and neutralinos/sneutrinos 
(d-f). 



in the previous section. This would amount to omitting factors of 1 in the brackets in the second Une 
of eq. (2.8) and in eq. (2.9) and to similar changes in the tadpole contributions. These changes do not 
affect, however, the interesting part of the threshold corrections which is not proportional to 5^^ . 

In the MSSM eq. (|1|) is replaced by -5"^^ [Tuo cos a /Ml + T^o sina/M|r] / Vu ■ This can be brought 
to a more convenient form by using the tree-level relations: cos^a/A/^ + sin^ Q/M|f = (sin^/3/A/§ + 



cos^ (i/M\)/ cos^ 2/3 and sin 



a cos a 



[l/Ml - l/Mjj) = -(sin/3cos/3/cos2 2/3)(l/Af| + l/M\). 
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Feynman diagrams describing contributions of chargino/charged slepton and neu- 
tralino/sneutrino sectors are shown in fig. ^. They give: 



7-charg 



3+Bk Y2j 



X 



{ml-Ml)B',iC,,E^) + BoiC,,E, 



(3.6) 



and 



rneutr 

'as 



92 + 9y ^AJ ^BJ* 



swZ]^ — cwZ'^J 



[m 



Ml)B',{N^,uj) + Bo{N,,uj] 



gl + glZ^'Z^-^*Z'^ (swZ]i - cwZ^ij) m^M^v ^j) 



(3.7) 



To check that the InQ dependence of the correction to the neutrino mass matrix in 
the MSSM matches the one following from the RGE it is more convenient to use the 
second approach described in the previous section and to assume that in the tree-level 
neutrino mass matrix is determined from the full 1-loop effective potential. Using then 



its RGE |5 



d 
di 



(■^Dl-loop - (f^)l-loop 



3 1 



{3.t 



and the MSSM RGE for C^^ (|r|) one finds 



a-loop^AB 



4M 



AM 



h~\oop_^ AB 



AM 



9 2 I 3 2 A 1 Q 

-292 + -^9y)ln- 



+ 



V 



1— loop 

4M 



(3.9) 



It is then easy to see that in the gauge ^vi/ = = (in which the effective potential and 
hence also t'l-ioop is defined) the dependence on InQ in the sum of corrections (PT^)-( PTT(]| ) 
and ( |3.4| )- ([3?^ ) cancels out with the InQ dependence in eq. ( |3.9| ). 



The final formula for the factor in the MSSM takes the form 



+ \ {g2Zf*Z'j* + ye,Zr-^'*Z'j*) {g2ZE'z'j + y.^Zl^^'Z 



1 + - n — ^ 

2 M^,^ 



X 
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_l_ 92 + 5'y ^AJ ^BJ 
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In 



+ terms proportional to 5"^^ 
where the functions / and g are 

1 a a? a 
f{a, b} = -- + + — In - 



+ g{mi,^,M^^) 

(3.10) 



-1 



a J a 
6 — a 6 



and satisfy f{a,a) = g{a,a) = 0. The terms proportional to the unit matrix are not 
interesting as they change only the overall scale of the neutrino masses and do not influence 
the mixing angles. 

Consider now the simplest limit Mh± = Mr± = = rric = mj\f. = Ms- Using the 
solution ( |1.4|) and writing for Q ^ Mz the factors as 

Ms' 



exp 



^In — 

167r2 Q 



(3.11) 



where tQ = (l/ievr^) ln(Mp/Q), tc, = (l/ievr^) ^(M^/M^) and i/^^ = igi/2)imljM^){l+ 
tan^/?), we have (up to the overall normalization) 

.AB,, / , \AB 



(mr) iQ) oc (mr) (Ms) 



1 



167r2^^ Q 



s 



Adding the 1-loop correction in the form 



lQ^2jAB ^^AbQ'^Is 



2 



l + tan^ p) 



+ 



+ OUMn2^ 



(3.12) 



1 2 -^S 

+ -0 y^^ In — — + terms proportional to o 

4t ^/ 



AB 



(3.13) 



we recover, as far the logarithms are concerned, the same result that is obtained with the 
running using the MSSM RGE from the scale Mp down to Msusy, and then the SM RGE 
to run down to the M^y scale [^]. There is however a nontrivial extra non- logarithmic 
piece 

,2 ^2 



2 M^\8 8 
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which is missed by the usual procedure. This piece is usually less important than the 
effects of the evolution from the scale Mp down to Msusy but for large tan (3 it is more 
important than the running from Msusy to the Mw scale with the SM RGE. 



4 Numerical analysis 

In the SM the most important correction which changes the matrix structure of the 
neutrino mass matrix can be incorporated by substituting 



where 



rth 



9 2 



167r2 2 



^-^ln^ + ^(..ln..) 



(4.1) 



(4.2) 



where xa = ^l^^/Myy. Since are proportional to the Yukawa couplings y'^^, this cor- 
rection cannot change qualitatively the results obtained by integrating the RGE and can 
be most easily taken into account by stopping the RG evolution of the Wilson coefficient 
C^^ at the scale Q = Mwe'^^'^^. The remaining corrections affect only the overall 
scale of the neutrino masses and therefore are not interesting in view of the unspecified 
magnitude of the mass M in eq. (|1.2| ). 

/M, 



In the MSSM the contribution of and to is also proportional to 5 



AB 



m 



W 



and cannot change qualitatively the results of the RG evolution.^ However the effects of 
the genuinely supersymmetric contribution I^^^ = I'^b^ + ^^b"^ can be important 

/My^. It has been 



AB 



m 



because unlike the SM case, it is not necessarily proportional to 5 
demonstrated ^ |T0[, that Pab ^ if large, can lead to relation between the mixing angles 
different than the one obtained at the infrared fixed point of the RGE [^]. The numerical 
estimates made in refs. |]^, ^, |TD[ relied however on approximating the corrections by 
the pure wino contribution to Ey'^(O). Here we analyze the dependence of on the 



MSSM parameters using the full expression ( p. 101) . For simplicity we will assume that the 
mixing of the left- and right-handed charged sleptons is negligible. 

We begin by considering fiavour conserving slepton mass matrices. In this case = 
^ABjsuhy i^gj^g^ygg ^i^Q matrices and Z^"^ are diagonal in the generation space. With 
no mixing of the left- and right-handed charged sleptons, the chargino and neutralino 
contribution can be simplified to 



1 2 
4^2 



In 



' Ml± 



^The same result can be always obtained from the running by slightly changing the value of tan 13. 
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Figure 6: Correction Pj^^^{Me) as a function of the left-handed charged slepton mass 
for chargino mass 150 (solid), 250 (dashed), 500 (dotted) and 800 (dot-dashed lines) for 
tan/3 = 2 and r = M2//X = +5 and —1. In lower panels we show the results of retaining 
only the wino (left) and wino and bino (right) contributions. 
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+ 



9i + 9Y 



— g2Z^Zj'mc, 



11.. ' V 



1 2 



z 



2i 



(4.3) 



where Mp± 



and M^± are the masses of the A-th. generation left- and right-handed 



charged sleptons, respectively. The contribution of the right-handed charged sleptons 
^j^usy-j^ (last line of i W^ ) is again proportional to the corresponding Yukawa coupling 
ul^. Hence, it too only slightly changes the effects of the RG evolution and can be 
neglected here. The remaining part (Z™**'*')^ of ([4.3|) depends on the mass of the charged 
slepton (the mass of the sneutrino is related to it by the underlying SUl{'2) symmetry: 



^LA 



^ +COS 2/?M^), tan P and the parameters of the chargino/neutralino sector: fi 
and M2 (as is customary, in the neutralino sector we take Mi ^ O.5M2). Figure |^ shows the 

E iriM^, 



results of the numerical evaluation of {I^ 



susy \ 



and compares it to the results 



of approximating P'^'^^{M^± ) by the pure (charged wino) or W^, and B (bino) 
contribution. The striking difference between the complete and approximate calculation is 
mainly due to the contribution of diagrams §3, c and e, f which give a negative contribution 
(third and fourth lines in eq. (|4.3| )) but are missed in the approximation. Although 
the absolute magnitude of the the correction P^^^{M^±^) does depend on tan/5 and the 
chargino composition, the differences P'^^^{M^± ) — P^^^{M^± ) (which are relevant for 
the changes of the neutrino mass matrix structure) are much less dependent on tan/5. 
They are however sensitive to the chargino (and neutralino) composition as is clear from 
the comparison of the two upper panels of Fig. 

In ref. [|] it has been observed, that if] |Jf > and J;;'^"^ - (/^""^ - P/) ^ 0, 

the masses of three neutrinos equal at the scale Mp can be split in agreement with the 
experimental information, provided the solar mixing angle is (very close to) maximal. 
For this mechanism to work — /™^^| ~ 10~^ is required. We can now improve 



the estimates made in ref. 



on the basis of the wino approximation. From Fig. ^ it 

ih by sub- 



^The renormalization group corrections (1.4), (1.5) can be incorporated in the formula 
stituting /' 



rth 



7-th 



I'^. One then has 



C^^(g) « C^^'iMp) + [P^iQ) - /^(Mf/Q)] C^^(Q) + C^'^iQ) [I'l^iQ) - IbHMf/Q)] . (4.4) 
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is clear that for M2I y^. ~ — 1 and lighter chargino mass 



1.2M^^) is sufficient to obtain r^^^ 



150 GeV the mass splitting 

~ 10-3 ( 10-3). 

For heavier charginos and/or M2//X positive obtaining — /^"^•*'| ~ 10^^ 
large mass splitting (or is even impossible to achieve). 



requires very 



X 10~^ 




X 10~ 


0.05 




0.05 












-0.05 




-0.05 


-0.1 




-0.1 


-0.15 




-0.15 


-0.2 




-0.2 


-0,25 




-0.25 


-0.3 




-0.3 


-0.35 




-0.35 


-0.4 




-0.4 


-0.45 


L r= + 5 J 


-0.45 


-0.5 




-0.5 




200 400 600 800 

(GeV) 



1000 



200 400 600 800 

M,. (GeV) 



1000 



Figure 7: Coefficient of the mass insertion as a function of the left-handed charged 
slepton mass for chargino mass 150 (solid), 250 (dashed), 500 (dotted) and 800 (dot- 
dashed lines) for tan/5 = 2 and r = M^j [i = —1 and +5. 



From the formula ( [4. 3] ) we can also quantify the magnitude of the off-diagonal cor- 
rections induced by the flavour mixing in the slepton mass matrices. Assuming that 
left-handed charged slepton masses are all approximately equal (which implies that sneu- 
trino masses are also all approximately equal), this is most easily done in the so-called 
mass insertion approximation ||16[|. In eq. ( [4.3|) terms involving e.g. charged left-handed 
sleptons can be written as 

ZfZ^''*H{Ml±) = ZfZ^''*H{Ml + {Ml± - Ml± )) 



k k 



6 



AB 



H{My - Ml^H'{Ml^J\ + Zt'Ml^Zi'*H'{Ml, 



{ML) 



AB 



H'{Ml± ) + terms proportional to (4.5) 



where Mav is the common mass of the left-handed charged sleptons and H is some function. 
We have also used the defining property of the matrices Zl. The mass insertion are defined 
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as the ratios of the off-diagonal elements (M|^)^^ of the mass squared matrix of sleptons 
to Ml^. In the case of flavour mixing only in the left-handed slepton sector[| we get for 



7-susy 



Tsusy 



iAB 



+ 



167r2 U 



9l 



2 



92 + 9yZn {swZjj - cwZfj^ mN^G{mlf^,M? 



(4.6) 



where 



F{a,b) 
G{a,b) 



6^ - 3ab 2a% 



iP 



iP 
ab 



a 



1 " 
3 ^"6 



1 ^ 



(4.7) 



We have used the fact that because of the underlying SUl(2) symmetry, mass insertions 
in the left-handed charged slepton sector and in the sneutrino sector are the same (and 
Ml± and M?^^ are related). In fig. |^ we plot the coefficient of as a function of M'i,± 
for several values of the chargino masses for tan/3 = 2. We see that for a fixed the 
biggest values of are obtained for M2/ ^ —1 and for rather large slepton to chargino 
mass ratio. In principle the mass insertion approximation should fail for \5xy\ ~ 0.1. In 
practice it works as an order of magnitude estimate even for \5xy\ ~ (^^^ error is then 
of order 25%). More accurate results can be always obtained from the general formula 



5 Conclusions 

We have computed the low energy threshold corrections to neutrino masses and mixings 
in the SM and in the MSSM. We have explicitly demonstrated that they are gauge in- 
dependent and stabilize the results with respect to the variation of the scale Q to which 
the relevant RGE is integrated from the high energy scale of the see-saw mechanism, thus 



clarifying the points raised in ref. [|TT 



^Flavour mixing in the right-handed slepton sector gives suppressed by UeAlJeB- Flavour non- 

diagonal entries in the soft supersymmetry breaking terms mixing charged left- and right-handed sleptons 
give proportional to 1722/63 and, hence, substantial only for large tan (3 values, i.e. only when the 

renormalization group corrections are dominant. 
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The general formulae for the corrections derived in this paper can be applied to 
various models predicting the neutrino masses and mixing. They can be used to quantify 
the slepton mass splitting and/or the amount of flavour violation in the slepton sector 
necessary to realize the specific mechanisms, investigated in papers ||, |lO|, allowing 
to obtain correct mass squared differences and mixing angles from initially equal neu- 
trino masses. They can find particularly interesting application in concrete models [0 
relating the see-saw mechanism generating neutrino masses to flavour non-conservation 
in the slepton sector. Finally, they will be indispensable for future precision tests of any 
quantitative theory of neutrino masses. 
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